SELF-AVERAGED SCALING LIMITS FOR RANDOM PARABOLIC WAVES 
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Abstract. We consider 6 types of scaling limits for the Wigner-Moyal equation of the parabolic 
waves in random media, the limiting cases of which include the radiative transfer limit, the diffusion 
limit and the white-noise limit. We show under fairly general assumptions on the random refractive 
index field that sufficient amount of medium diversity (thus excluding the white-noise limit) leads 
to statistical stability or self-averaging in the sense that the limiting law is deterministic and is 
governed by various transport equations depending on the specific scaling involved. We obtain 6 
different radiative transfer equations as limits. 
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1. Introduction 
The celebrated Schrodinger equation 



d^> ti 2 

ih— + —AV + (rV(t,x)y = 0, *(0,x) = * o (x) 
O ■ ot 2m 



describes the evolution of the wave function ^ of a quantum spin-less particle in a potential aV 
where a is the typical size of the variation. 

A similar equation called the parabolic wave equation is also widely used to describe the propa- 
gation of the modulation of low intensity wave beam in turbulent or turbid media in the forward 
scattering approximation of the full wave equation 19 j. In this connection the refractive index 
fluctuation plays the role of the potential in the equation. Nondimensionalized with respect to the 
propagation distances in the longitudinal and transverse directions, L z and L x , respectively, the 
parabolic wave equation for the modulation function reads 



(1) ik- 1 L~ l ^- + 2- 1 k- 2 L- 2 ^ + aV(zL z ,*L x )V = 0, *(0,x) = * (x) 



where k is the carrier wavenumber, ^ the amplitude modulation and A the Laplacian operator in 
the transverse coordinates x. Here we have assumed the random media has a constant background. 
In the sequel we will adopt the notation of (^Q). 

In this paper we study the scaling regimes where the wave field experiences both longitudinal 
and transverse diversity of the random medium, represented by V, whose fluctuation is assumed to 
be weak. This gives rise to random spread of wave energy in the transverse directions. To fix the 
idea, let us choose the units of the longitudinal and transverse coordinates such that the spectral 
support of the function V is O(l) while 

(2) L z ~ L x » 1, 

namely L z and L x are much larger than the correlation lengths of the medium in the longitudinal 
and transverse directions, respectively. There is no loss of generality in the choice of the hyperbolic 
scaling ((2j) (cf. Remark 2 below). Depending on the actual length scales and anisotropy of the 
medium, we adjust the intensity a of the medium fluctuations to obtain a nontrivial limit. 
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2. WlGNER DISTRIBUTION AND TIME REVERSAL 

There has been a surge of interest in the radiative transfer limit in terms of the Wigner distribu- 
tion (see below) because of its application to the spectacular phenomena related to time-reversal 
(or phase-conjugate) mirrors [I], [3], [El, [TB] - 

The Wigner transform or distribution of the wave function \E f is defined as 



W(z,x,p) = *J e ^™*(z,xi+Z)**(z,x-Z)dy 



1 / /' e -iP-(yi-y2) ( 5( x _yi+Z2) p ( yi)y2 ) dyidy2 



(2vr) c 

where p(yi,y2) = ^ r (yi)^ r *(y2) is called the two-point function or the density matrix. As appar- 
ent from the definition the Wigner function contains all the information about p. The Wigner 
distribution have the simple properties 

J W(z,x,p)dxdp = ||*|||, || W || 2 = (2?r)- d / 2 

This is the case of pure state Wigner distribution. What is the more pertinent for us is the so 
called mixed state Wigner distribution. 

Let us briefly review how a mixed state Wigner distribution arises in the time-reversal operation. 
Let Gff(0, x, z, y) be the Green's function, with the point source located at (z,y), for the reduced 
wave (Helmholtz) equation for which the Schrodinger equation is an approximation. By the self- 
adjointness of the Helmholtz equation, Gh satisfies the symmetry property 

G H (0,x,z,y) = G H (z,y,Q,x). 

The wave field ^ m received at the mirror is given by 

) / Gtf(0,X m ,2!,X s )*o(Xs)cix s 

= XA(x m ) J G H (z,x s ,0,x m )^ (x s )dx s 

where \A is the aperture function of the phase-conjugating mirror A. 

After phase conjugation and back-propagation we have at the source plane the wave field 

fy B (z,x;k) = j G H (z,x,0,x m )G H {z,x s ,0,x m )xA(xm)^o(xs)dxmdx s . 

In the parabolic approximations the Green's function Gh{z, x, 0, y) is approximated by e lkz Gs{z, x, y) 
where Gs(z, x, y) is the propagator of the Schrodinger equation. Making the approximation in the 
above expression for the back-propagated field we obtain 

$> B (z,x;k) = J Gs(z,x,x m )Gs(z,x s ,x m )^ (xs)XA(xm)dx m dx s 

(3) = J e^-^W(z, p)^(^)dpdx s 

where the Wigner distribution W is given by 

(4) W(z,x,p) 

= J e~ ip ' y Gs(z,x + y/2,x m )G s (z,x - y/2,x m )xA(x m )dydx m . 

This is a mixed-state Wigner distribution. In general, the integral in © should be interpreted in 
the distributional sense. 



The Wigner distribution in (J3J) has the initial condition W 
(5) W(0,x,p) = ^0 

and can be treated as a generalized function on M. 2d . Indeed, for any O £ C£°(R d ) we have 



(6) (q> B ,Q) = J J W(z,r,p)e(r,p)drdp 



where the function is defined as 



0(r, p) = 2 d J e(y)e i2p -(y- r )/T^ (2r - y)dy. 



If for instance G C£°(K ) then it is easy to see @(y, p) is compactly supported in y £ R" 
and decays rapidly (faster than any power) in p € M. d . As a result we can always approximate to 
arbitrary accuracy the distributional initial data such as © by square-integrable initial data. 

The fluctuations of the back-propagated wave field is thus determined by the fluctuations of 
the Wigner distribution. The statistical stability or self-averaging of the Wigner distribution in 
turn explains, modulo the scaling limit, the persistence and stability of the super-focusing of the 
time-reversed, back-propagated wave field observed experimentally and numerically. 

Our main results show that under various scaling limits, sufficient amount of spatial-transverse 
diversity experienced by the propagating wave pulse results in self-averaging and deterministic 
limiting laws. 

From the perspective of the quantum stochastic dynamics in a random environment, our results 
say that due to the spatio- temporal diversity experienced by the wave function of the quantum 
particle the quantum dynamics has in the scaling limit a classical probabilistic description which 
is independent of the particular realization of the environment. The transition from a unitary 
evolution to an irreversible process is of course the outcome of the phase-space coarse-graining by 
the test functions. The results presented below are a rigorous demonstration of decoherence, a 
mechanism believed to be responsible for the emergence of the classical world from the quantum 
one P2], H3- 

3. Assumptions 

We assume V z (x) = V(z, x) is x-homogeneous square-integrable process with the (partial) spec- 
tral measure V(z,dq). which is an orthogonal random measure 

Eft(z,dp)V(z,dq)] =5(p + q)$ (p) d P dc l 
and gives rise to the (partial) spectral representation of the refractive index field 

V z (x) = V(z,x) = J exp (ip • x)V(z, dp). 

In case that V(x),x € is x-homogeneous square-integrable random field with the full spectral 

density given by 3>(£,k) we have the following relation 

$o(p) = J Hw,p)dw. 
We also have the following relation between the partial and full spectral measures 

V z (dp) = J e izw V{dw,dp) 
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such that 



E[y 2 (dp)F s (dq)] 



j e^ s ~ z ^ w ^(w,p) dw 5(p + q) dpdq 
<5(s — z, p)5(p + q) dp dq 



where 



e isw $(u;,p) dw. 



Since $(k) 
(7) 



Vio, p 



$(s,p) = 

<&(— k),Vk £ we may assume in this case that 

<&(u;, q) = <f>(-w, p) = -p) = $(-w, -p), 

so that <l(s,p) is real-valued and <£(s,p) = $(— s,p). 

For simplicity of the analysis we assume that the spectrum is smooth and has a compact 
support. We also assume the sample field is almost surely smooth and its spatial derivatives all 
have finite moments. The conditions of smoothness and having a compact support are unnecessary; 
they are assumed here in order to make simple the key estimate (Proposition 5, 6, 7, 8). 

Let T z and be the sigma-algebras generated by {V s : Vs < z} and {V s : Vs > z}, respectively. 
Define the correlation coefficient 



(8) 



p{t) = sup 

E[h]=0,E[h 2 l = l z +* 

E[ 9 ]=0,E[g2] = l 



sup 



E [hg] . 



Assumption 1. The correlation coefficient p(t) is integrable 

When V z is a Gaussian process, the correlation coefficient p(t) equals the linear correlation 
coefficient r(t) which has the following useful expression 



(9) 



r{t) = sup I &(t-Ti-T 2 ,k)gi(Ti,k)g 2 (T2,'k)dkdTidT2 

91,92 . 



where the supremum is taken over all g±,g 2 £ L 2 
satisfy the constraint 



_1 ) which are supported on (— oo,0] x R'' and 



(10) j $(t-t',k)g 1 (t,k)g 1 (t',k)dtdt'dk = I 4>(t - t', k)g 2 (t, k)g 2 (t' , k)dtdt'dk = 1. 

There are various criteria for the decay rate of the linear correlation coefficients in the literature. For 
example, according to Chapter V, Theorem 8 (after adaptation to the t-continuous version), 
even under the restrictive condition that $ € there is a large class of Gaussian processes whose 
correlation coefficient decays faster than any power law. 

Secondly, we assume a 6-th order quasi-Gaussian property: Let 

C/ s 1 (x) = T/ s (x), C/ s 2 (x)=E,[y s ](x), s>z. 

Assumption 2. For any choices of o~j G {l,2},j = 1,2, ...,N and a set of linear operators {Tj}, 
there exists a finite constant C 



E 



E 



A' 



ll^rl'Mx,) 

3=1 



0, iV = 3,5 



N 



He [r m [/^(x m )r n [/f : (x n )] 



N = 4,6 



where the summation is over all possible pairings {mn} among {1,2, ...,N}. 
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Finally we assume 

Assumption 3. There exists a constant C such that for Theorem 1, 2, 3 (i), (Hi) and 4 (i), (Hi) 

(11) limE[sup \\C £ z e\\l] < -E||£^||l, G C C °°(M M ), W < oo; 

e^O Z<ZQ e 

and for Theorem 3 (ii) and 4 (ii) 

limE[supl|£!e|||] < — E||£!0||1, V0 G C c °°(M 2d ), Vz < oo; 
z<z £ a 



c 

e^o L Z<ZQ " " - ~ e 



limE[sup \\C £ z C £ z e\\l] < -^E\\C £ z C e z 0\\2, V# G C~(R 2d ), Vz < oo; 



where C e z is defined, respectively, by \6Q(l . h92\). H0(]l\) and h!21\) and a G (0, 1) as specified in the 
statements of the theorems. 

Assumption 3 is readily satisfied for Gaussian random fields. This can be seen by first observing 
that C £ z 6 is a Gaussian process and C £ z C e z 9 is a x 2 -process and, secondly, by an application of BorelPs 
inequality that the supremum over z < zq inside the expectation can be over-estimated by a 
log (1/e) factor for excursion on the scale of any power of 1/e: 

(12) E[sup \\£ £ z ef 2 ] < Clog(-)E||£^|| 2 ; 

z<zo £ 

(13) E[sup \\£ z £ie\\l} < Clog 2 (-)E||r 2 £^||i. 

Z<ZQ £ 

4. Main results 

In the standard scaling, we set 

(14) L z = L x = » 1, a = e. 

To describe the small scale wave energy we consider the scaled version of the Wigner distribution 
W £ (z,x,p) = e -^(z^ + £ ^W(z^- £ ^)dy 

1 f f ,-ip-( yi -y 2 )/ £ 2 r/„ yi+y2- 



(15) = (2^F 7 J e ~ iP ' m ' y2)/e 8(x-^^)p(y l ,y 2 )dy 1 dy 2 . 

The Wigner distribution W £ has a limit as certain measure, the Wigner measure, introduced in |15j . 
But as remarked in the introduction, we always consider uniformly L 2 initial condition induced by 
a mixed-state density matrix p. 

The Wigner distribution satisfies the Wigner- Moyal equation 

(16) + I • VW £ + -C £ Z W £ Z = 

with W z (x, p) = W £ (z,x, p). Here the integral operator C £ is given by 

(17) £ £ z W £ (x, x, p) = if e** [W|(x, p + q/2) - W £ z (x, p - q/2)] V(^, dq), x = xe" 2 " 
with a = 1. 

The more general case with a G (0, 1) can be derived from a somewhat different scaling (cf. 
the scaling leading to Theorem 2): We probe a highly anisotropic medium V(z, e 2 ~ 2a x) with the 
strength 



with a wave beam composed of waves of lengths comparable to that of the medium, so we replace 
k by ke 2 ~ 2a : 

(18) k — » ke 2 - 2a 

in the parabolic wave equation. We then use the following definition of the Wigner distribution to 
resolve the wave energy: 

(19) W £ (z, X ,p) = 1 le-^f^x + ^z.x-^y. 

The difference in scaling between (fT$|) and (fTH|) is, of course, due to the rescaling of coordinates 

HUD; 

Since the proof of convergence is the same for a € (0, 1] they are treated together. Eq. (|2T)|) 
and its variants studied in the sequel are understood in the weak sense and we consider their 
weak solutions with the test function space C^°(R 2d ): To find W £ G L 2 ([0, oo); L 2 (R 2d )) such that 
||W z e || 2 < ||Wo|| 2 ,V^ > 0, and 

(20) (W £ z , 9) - (Wo, e) = k- 1 f (W £ s , p • V x #) ds + - f (W £ ,£ £ s 6) ds. 

Jo £ Jo 

We shall use the notation 

V £ (dq L )=V(^,dq), y/(x)=F(4,x). 
Taking the partial inverse Fourier transform 

F^e^y)^ J e^(x,p)dp 

we see that J-^ 1 C £ 8(x,y) acts in the following completely local manner 

(21) ^/^(x.x.y) = -i<5 e F/(x,y)^ 2 -^(x,y) 
where 

(22) 5 £ V z £ (it, y) = F/(x + y/2) - F/(x - y/2). 

The operator £| is skew-symmetric and real (i.e. mapping real-valued functions to real-valued 
functions). 

Since our results do not depend on the transverse dimension d we hereafter take it to be any 
positive integer. 

Remark 1. Since Eq. \20]) is linear , the existence of weak solutions can be established straightfor- 
wardly by the weak-* compactness argument. Let us briefly comment on this. First, we introduce 
truncation N < oo 

V N {z,x) = V(z,x), \V(z, X )\<N 

and zero otherwise. Clearly, for such bounded Vn the corresponding operator C £ z is a bounded 
self-adjoint operator on L 2 (R 2d ). Hence the corresponding Wigner-Moyal equation preserves the 
L 2 -norm of the initial data and produces a sequence of L 2 -bounded weak solutions. Passing to the 
limit N — ► oo we obtain a L 2 -weak solution for the original Wigner-Moyal equation ifV is locally 
square-integrable as is assumed here. However, due to the weak limiting procedure, there is no 
guarantee that the l?-norm of the initial data is preserved in the limit. 

We will not address the uniqueness of solution for the Wigner-Moyal equation \2U}) but we will 
show that as e — > any sequence of weak solutions to eq. \2U]) converges in a suitable sense to the 
unique solution of a deterministic transport equation. 

We state our first result in the following theorem. 
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Theorem 1. Let <3? G and Assumption 1, 2, 3 be satisfied. Then the weak solution of the 
Wigner-Moyal equation \2U\) . Ji7| ) with the initial condition Wq G L 2 (M? d ) converges in probability 
as the distribution-valued process to the deterministic limit given by the weak solution W z of the 
radiative transfer equation 

+ f •V^(x,p)=^(x,p) 

with the initial condition Wo and one of the following scattering operator C: 
Case (i): < a < 1, 

(23) £W z (x,p)=2ir J <D(0,q-p)[V^(x,q)-^(x,p)]dq; 
Case (ii): a = I, 

(24) £W 2 (x,p) = 27r J $( lq|2 ~ |P|2 , q - p)[W 2 (x, q) - W z (x, p)]rfq; 

Case (iii): a > 1, 

(25) CW Z = 0. 

The case of a = corresponds to the so-called white-noise scaling whose limit is a Markovian 
process Eq. (|24jl has recently been obtained in [3] for strongly mixing z-Markovian refractive 
index fields with a bounded generator. 

In order to obtain a nontrivial scattering kernel for a > 1 we need to boost up the intensity of 
V (cf. Theorem 3). 

Next we consider a second type of scaling limits which starts with the highly anisotropic medium 
V(z,e 2 - 2a x). We then set 

(26) L x = L z = e~ 2 , a = e 2a ~ 1 , 0<a<l 
under which the parabolic wave equation becomes 

(27) ik- 1 ^- + 2- 1 k~ 2 e 2 A* £ + e 2a ~ 3 V(ze- 2 ,xe- 2a )y £ = 0, ^ £ (0, x) = ¥„(x). 

The radiative transfer scaling (|14|) is the limiting case a = 1. The time-evolution of the Wigner 
function (|15j) is governed by the Wigner-Moyal equation (|2()|) with the following operator C e z 

(28) £^(x,x,p) 

= i j e iq 'V a ~ 2 [W z £ (x, p + £ 2 ~ 2a q/2) - W|(x, p - £ 2 ~ 2o q/2)] F-f(dq), x = xe" 2 ". 

The partial Fourier transform of is now given by (|21j) with the following (5 £ T4 

(29) 6 £ V*(Z, y) = 6 2 "" 2 [y/(x + ye 2 " 2 "^) - K e (x - y £ 2 - 2 72)] . 
We now state the result for the scaling limit (|26j). (|27|). 

Theorem 2. Let < a < 1. Let $ G C£° and Assumption 1, 2, 3 be satisfied. Then the weak 
solution Wg of the Wigner-Moyal equation K2U\) . 128\) with the initial condition Wq G L 2 (M. 2d ) 
converges in probability as the distribution-valued process to the deterministic limit given by the 
weak solution W z of the following advection- diffusion equations with the initial condition Wo: 

(30) ^ + R.yw z = k 2 V p -~DW p W z , 

with one of the following diffusion tensors D. 



• a G (0,1): 

(31) D = tt J $(0,q)q® qdq; 

• a > 1: 

D = 0. 

For a = 1 the limit is the same as that in Theorem 1 Case (ii). a = gives rise to the white- noise 
limit for the Liouville equation. The advection-diffusion equation (J30J) can be obtained from (|23|) 
under the diffusion limit of the latter. 

Let us consider yet another type of scaling limit parametrized by (3. We first assume a highly 
anisotropic medium V{e 2 ~ 2 P z,~x) and set 

(32) L x = L z = e~ 2 , a = e 

i.e. the radiative transfer scaling. The Schrodinger equation then becomes 
r)\Er e 

(33) ik- 1 —- + 2- 1 k- 2 e 2 ^ £ + e- l V(z£- 2p ^e- 2 )^ £ = 0, f £ (0, x) = * Q (x), (3 < 1 

the corresponding Wigner-Moyal equation is 

f)W £ n h 

(34) + P . VW e + = 

with 

(35) £fT^(x,x,p) 

= ^-^e^tW^^p + q/^-W-^^p-q/^y^^q), x = x^ 2 , (3 < 1. 

Eq. 1)351 is a borderline case of the following family of scaling limits. Let us consider probing 
an anisotropic medium V{e 2 ~ 2 ^ 'z, £ 2_2q x), a > 0, with a wave beam composed of waves of lengths 
comparable to that of the medium, so we switch to ()18|) and (|19|) for the formulation of scaling 
limits. 

Two situations arise: Case (i) a < j3 and Case (ii) a > (3. 

In the first case a < j3 we set the strength of the medium fluctuation to be 

a = e 2a -P. 

The resulting equation is ()34|) with 

(36) £f^(x,x, P ) = iy"e^*[W'|(x,p + q/2)-W^(x,p-q/2)]F( i ^,dq), x = xe" 2 " 
In the second case a > (3 we set the strength of the medium fluctuation to be 



<t = e . 



The resulting equation is ()34|) with 
(37) £|W£(x,x,p) 

= ^- a y*e^[W|(x )P + q /2)-H^(x,p- q /2)]T/(-^,d q ), x = xe" 2 ". 

We have the following theorem. 

Theorem 3. Let a,f3 > 0. Let G and Assumption 1, 2, 3 be satisfied. Then the weak 
solution W £ of the Wigner-Moyal equation J^| ) with Case (i): or Case (ii): h'Jb)) and the 

initial condition Wo G L 2 (M. 2d ) converges in probability as the distribution-valued process to the 



deterministic limit given by the weak solution W z of the transport equation with the initial condition 
W : 

(38) ^ + R. VWz = k 2cW z , 

with one of the following the scattering operators C 
Case (i): a < (3 

(39) £W z (x,p) = 2tt f *(0,q-p)[W,(x,q)-Wi(x,p)] dq. 
Case (ii): 1 < or/0 < 4/3, d > 3 

(40) £W z (x,p)=2tt y ^( |q|2 ~ |p|2 ) y $(to,q-p)du; [W,(x, q) - W„(x, p)] dq. 
Case (iii): a = [3 

(41) £W,(x,p)=2tt y $( |q|2 ~ |P|2 , q - p) [W 2 (x, q) - W 2 (x, p)] dq. 

Theorem 3 (i) probably holds for d = 2 and a/ (3 > 4/3 but we do not pursue it here in order to 
keep the argument as simple as possible. 

Earlier JH]) |Ej have established the convergence of the mean field EW| for z-independent Gauss- 
ian media and d > 3. Their transport equation can be viewed as a limiting case of (]38[) in which 
$(<!;, k) is a 5-function concentrated at £ = 0. See also \17\ for mean-field results for z-finitely 
dependent potentials. 

Unlike the transport equations (|24"]) , p3|) , the scattering kernel l|4*0"|) is elastic in the sense that it 
preserves the kinetic energy of the scattered particle so that the incoming and outgoing momenta 
q, p have the same magnitude. 

Finally let us consider two other types of scaling limit starting with the slowly- varying refractive 
index field V(e 2 ~~ 2/3 z, e 2_2Q x), a, (3 € (0, 1). In the first case 

(42) (3 > a, < a < 1 
we set 

(43) L x = L z = e~\ a = e 2a ~ 13 
under which we have the following parabolic wave equation 

(44) ik- 1 ^- + 2- 1 k- 2 e 2 A^ £ + e 2a ^- V(ze" 2/3 , xe~ 2a )^ £ = 0, * £ (0,x) = * (x), 

and the corresponding Wigner-Moyal equation (|34|) with 
(45£|Wj(x,x,p) 

: y e *»V«-2 [WJ(x,p + e 2 - 2a q/2) - Wj(x,p - e 2 - 2Q q/2)] V(^,d q ), x = xe" 2 ". 

In the second case 

(46) a > fi, < a < 1, 
we set 

(47) L x = L z = e- 2 , a = e a . 
After rescaling the parabolic wave equation reads 

(48) ik' 1 ^- + 2- 1 k~ 2 e 2 Ay £ + e a - 2 V(ze- 2f3 ,xe- 2a )y £ = 0, * £ (0, x) = f (x), 
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and the corresponding Wigner-Moyal equation takes the form of (]34)l with 
(4££gW|(x,x,p) 

= ie p - a J e iq 'V Q - 2 [Wf(x,p + e 2 - 2Q q/2) - W|(x, p - e 2 " 2a q/2)] F(-^, dq), x = xe" 2 " 

Theorem 4. Xei a, (3 € (0,1). Let <& € and Assumption 1, 2, 3 be satisfied. Then the weak 
solution W% of the Wigner-Moyal equation \34\ ) with Case (i): \4-5\j or Case (ii): and the 

initial condition Wo G L 2 (M. 2d ) converges in probability as the distribution-valued process to the 
deterministic limit given by the weak solution W z of the advection- diffusion equation \3U\) with the 
following diffusion tensors: 

Case m- (j42]L (j43j): 

(50) D = tt j $(0,q)q®qdq. 

Case (ii)-d46j), (JUJ): d > 3, 1< a/(3 < 4/3, 



(51) D(p) = irk\p\- 1 J j <S>(w,px)dw 



p_L <g> p±dpj 



where p± G R rf 1 , p^ • p = 0. 
Case (iii): a = (3 



dq 



(52) D(p) = tt J $(fc 1 p ■ q,q)q<g>q 

The advection-diffusion equation with (|50|). (|51|1 and 1)52(1 are the diffusion limit of the transport 
equations (|40|1 and 1)24)1 . respectively. The limiting case of a = gives rise to the white-noise model 
of the Liouville equation jH]. We believe that the result for Case (ii) can be extended to d = 2 and 
P/a G (0,1). 

Remark 2. Taken together, our results have roughly covered all the super-parabolic scaling 

L x » y/lT,. 

To see this, let us set 

L Z ~L% = e~ 27 , < 7 < 2 
and define the following Wigner transform 

~2,, ^2 



1 

(27T 



!_ /" /'.-p-(yi-y 2 )/ £ - 2 ( 5 (x _yi±^ )p(yi)y2)ciyi(iy2 



wrf/i i/ie new parameter 



and analyze analogously the preceding scaling limits as parametrized by e. For a systematic treat- 
ment of scaling limits resulting in a transport equation, see • 



Our approach is to use the conditional shift [13] to formulate the corresponding martingale 
problem parametrized by e and adapt the perturbed test function technique to the probabilistic 
setting to establish the convergence of the martingales. It then turns out that after subtracting the 
drift and the Stratonovich correction term the limiting martingale has null quadratic variation (see 
Proposition EJ) implying that the limit is deterministic. The perturbed test functions constructed 
here (see e.g. (|71|). (jSTJl and (|82)l) are related to those in [2], [3] but our analysis is carried out in a 
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more general framework as formulated in l B| and provides a unified treatment of a range of scaling 
limits from the radiative transfer to the diffusion limit and the white-noise limit. 



(53) 



5. Martingale formulation 
We consider the weak formulation of the Wigner-Moyal equation: 

-l 



[(WZ,6)-{W ,0)] 



k' 



f (W £ ,p-V9)ds + - [ {W £ s ,C £ z 6)ds 
o £ Jo 



for any test function 6 E C^°(M. 2d ). which is a dense subspace in L 2 (M. 2d ). The tightness result (see 
below) implies that for L 2 initial data the limiting measure P is supported in L 2 ([0, z }; L 2 (R 2d )). 

For tightness as well as identification of the limit, the following infinitesimal operator A £ will 
play an important role. Let V£ = V(z/e 2 ,-). Let T z be the a-algebras generated by {V £ , s < z] 
and W z the corresponding conditional expectation w.r.t. T z . Let M £ be the space of measurable 
function adapted to {JRj, Vi} such that sup 2<20 E|/(z)| < oo. We say /(•) € V(A £ ), the domain of 
A £ , and A £ f = gif f,g£M £ and for f s (z) = 5- l [W z f{z + 5)- f(z)] we have 



supE\f\z) 

z,S 

limE\f s (z)-g(z) 

o— >U 



< OO 

= 0, Vz. 



Consider a special class of admissible functions f(z) = <p((W £ ,6)), f'(z) = 4>'((W £ , 6)), V0 € C c 
we have the following expression from ()53|) and the chain rule 



(54) A £ f(z) = f'(z 
In case of the test function 8 that is also a functional of the media we have 

(55) A £ f(z) = f'(z 



\ (W £ , p-V9) + ^ (W £ z ,C £ z 0) + (W £ z , A £ 



and when 6 depends explicitly on the fast spatial variable 

x = x/e 2a 

the gradient V is conveniently decomposed into the gradient w.r.t. the slow variable V x and that 
w.r.t. the fast variable V-z 



V = V x + £~ 2Q Vx. 



A main property of A £ is that 



(56) 

Also, 

(57) 



/(*) 



[ A £ f{s)ds is a .Ff -martingale, V/ G V{A £ ). 
Jo 



E £ f(z)-f(s) = J\ £ A £ f(r)d- 



t Vs < z a.s. 



(see ^H])- We denote by A the infinitesimal operator corresponding to the unsealed process V z {-) 
V(z,-). 

li 



6. Proof of Theorem 1 
6.1. Tightness. In the sequel we will adopt the following notation 

/Oz) = 0(<w|,0>), f(z) = <t>\{w s z ,e)), f"( 2 ) = <t>"((wi,e)), V0ec c 

Namely, the prime stands for the differentiation w.r.t. the original argument (not t) of /, /' etc. 

Let -D([0, oo); L^,(M 2d )) be the L 2 -valued right continuous processes with left limits endowed with 
the Skorohod topology. A family of processes {W £ , < e < 1} C D([0, oo); Lf u (M 2d )) is tight if and 
only if the family of processes {{W £ ,9) ,0 < e < 1} C D([0,oo); Ll(R 2d )) is tight for all 9 G C c °° 
jlUj . We use the tightness criterion of fT3, (Chap. 3, Theorem 4), namely, we will prove: Firstly, 

(58) lim limsupP{sup | {W £ , 9) \ > N} = 0, Vz < oo. 

N^oo £ ^q z <z 

Secondly, for each <j) G C°°(R) there is a sequence f £ (z) G V(A £ ) such that for each z < oo 
{A £ f £ (z), 0<e<l, 0<z< zq} is uniformly integrable and 

(59) \im¥{sup \f £ (z) -<p{{W £ ,9})\> 5} = 0, V«5 > 0. 

e^O 2<20 

Then it follows that the laws of {{W £ , 9) , < e < 1} are tight in the space of D([0, oo); R). 
First, condition (|58[) is satisfied because the L 2 -norm is uniformly bounded. 
Let 



(60) ^(x,x,p) 



J e 4q ' i [0(x, p + q/2) - 0(x, p - ci/2)}e lk ~ 1{s ~ z) ^ /£2a W z V £ {dci)ds 



which has a compact support since both 9 and $ do. Note that the operator C £ z maps a real-valued 
function 9 to a real-valued function. 

Lemma 1. 









" /'OO 




(61) E 




2 (x,p) < 


/ p(s)ds 
Jo 
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(x, p + q/2) - 0(x, p - q/2)] 2 q)o?^q 

which has an e-uniformly bounded support and is bounded uniformly in x, p, e. 
Proof. Consider the following trial functions in the definition of the maximal correlation coefficient 
h = h s (x,p) 

= i j e iq - X£ ~ 2 >(x,p + q/2) - 0(x,p- q/2)]e i ^ 1(s - z)p - q/E2 "E^/(dq) 
9 = 5t(x,p) 

= i j e^^I^x^ + q^-^x^-q^le^ 1 ^^ 2 ^/^) 
It is easy to see that 

h s G L 2 (P,n,T £ -2 z ) 



gt G £L 2 (P,n,F 



e-n 



and their second moments are uniformly bounded in x, p, e since 

(62) E[/i 2 ](x,p) < E[ 5 2 ](x,p) 

(63) E[ 5s 2 ](x,p) = y"[fl(x,p + q/2)-e(x,p-q/2)] 2 $(C,qKdq 
which is uniformly bounded for any integrable spectral density <3?. 
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From the definition © we have 

|E[/i a (x,p)/i t (y,q)]| = \E[h.gt]\ < p(e~ 2 (t - z))® 1 / 2 [h 2 a (x, p)] E 1 / 2 [ 5 2 (y,q)] . 
Hence by setting s = f,x = y,p = q first and the Cauchy-Schwartz inequality we have 

E[/i 2 (x,p)] < p 2 (e- 2 (s-z))E[g 2 ( X ,p)] 
|E[^(x,p)^(y,q)]| < p(e- 2 (t-z))p( e - 2 ( S -z))E 1 /2[ 5 2( X) p)] E i/2[ 5 2 (yiq) ] ) Vs ,i>z,Vx,y. 
Hence 



noo 
E[/i,(x,p) ft (x,p)]cfed4 < E[ 5t 2 ](x,p) 



p(s)ds 



which together with (|6*5)) yields (|6*T]) . 
Corollary 1. 

(64) E 



□ 



P • V x £* 

< 



- 2 

p(s)ds 







y [p • V x 0(x, p + q/2) - p • V x 0(x, p - q/2)] 2 q)d£dq 

which has an e-uniformly bounded support and is bounded uniformly in x, p, e. 
Inequality Q64[) can be obtained from the expression 
p • V x £*0(x,x,p) 



IE 



J°° J e *»-*p • V x [0(x, p + q/2) - 0(x, p - q /2)y k -^ s ~ z ^ 2a ElV^d<i)d S 



as in Lemma n 

The main property of C e z 9 is that it solves the corrector equation 



(65) 



e~ 2a T- V x + ^ £ 
k 



CIO = e~ l L 



2 re. 



Eq. (|65|) can also be solved by using (|21jl. yielding the solution 

/oo _ 
e _ 4£ - 2 . fe -i (s _ z)Vy .v, [ E « (x,X,y)dfl 

where 

6 E V z %x, y) = K £ (x + y/2) - V/(x - y/2). 
Recall that V x and V x are the gradients w.r.t. the fast variable x and the slow variable x, 
respectively. 

We will need to estimate the iteration of C% and C% : 
£|£|0(x,x,p) 

(67) 

£^0(x,x,p) 
(68) 



y/(dq)E £ [l//(dq / )]e 4q - i e iq '- i e ifc " 1( ^^ p - q/£2Q 
{[0(x, p + q'/2 + q/2) - 0(x, p + q'/2 - q^^^-^q'-q/^ 2 "*) 
-[0(x, p - q'/2 + q/2) - 0(x, p - q'/2 - q /2)] e - ifc " 1 ( s -^'-q/( 2e2Q ) J ds 

I W z [Vf (dq)]E £ [V t £ (dq')] e iq - V q ' ' V fc ~ 1 ( s - z )p-^ 2q e ifc " 1 (t-*W/e a 
{ [0(x, p + q'/2 + q/2) - 0(x, p + q'/2 - g^)]^ 1 ^^'^ 2 ") 

-[0(x, p - q'/2 + q/2) - 0(x, p - q'/2 - q /2)] e - ifc " 1 («-^)q'-q/(2 E 2a ) J dsdt 
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Both have a compact support. Their second moments can be estimated as in Lemmanby using the 
6-th order quasi-Gaussian property (Assumption 2). In order to carry out the same argument we 
need to approximate the terms of non-product form such as 0(x, p + q'/2 + q/2)e ifc ( s_2: ) c i'' < i/( 2e a ) 
by sum of product of functions of variables that are statistically coupled in the quasi-Gaussian 
pairing. 

Since we do not need the pointwise estimate such as stated in Lemma|I]we shall demonstrate a 
simpler approach based on the inverse Fourier transform: 

(69)^ [Clltej (^y) = e- 2 J tf^e"*" " fc " (-*> v *- v * [E^V^O] (x, x, y)ds 



(70)^ 2 - 1 {£^}(x ! x,y) 



-ie- 2a k- 1 (t-z)V 



[E z [5 £ V £ ]T 2 l 6]} {^±,y)dsdt. 



Lemma 2. 



E\\£ e z £l9\\i < SCl^J^ p(s)dsj E[V Z )' J [0(x,p + q/2)-0(x )P -q/2)]^(^q)^xdqdp 

E\\C £ z C £ z e\\ 2 2 < 8C^J p(s)ds^j E[yj 2 y[0(x,p + q/2)-0(x,p-q/2)] 2 $(e !q Kdxdqdp 
for some constant C independent of e. 

Proof. The calculation for C £ Z C £ Z 9 is simpler. Let us consider C e z C e z 9. 
By the Parseval theorem and the unitarity of exp (irV y ■ V x ), t G M, 

E\\t z cie\\l = E\\T 2 l t z cie\\l 

= e~ 8 J E j^V^ 2 ^ 1 ^^-^ [E z [5 £ V t £ ]e- l£ ' 2ak ~^ s ~ z ^^- [E^V^T^O] } (x, x, y)dsdt 



y^{E z [5 £ V t £ 



e] ie- 2a k- 1 (s'-z)VyS/ i 



[E z [8 e V £ ,]^ x e] } (x, x, y)ds'dt'\ d*dy 



E 



ie- 2a k- 1 (t-t')/2\7 : 



' v * [E z [5 £ V £ ]e- ie - aa *- 1 (-*)VyV* [e z [5 e V*\F 2 x 6]) (x,St,y)dsdt 



ie- 2o 'k- 1 (t-t')/2\7 



y^^E^V^- 2 ^- 1 ^-^*^ [ Ez [6 £ V s ^ 2 -^}(x,x,y)ds'dAdxdy 



e- 8 I E 



[E z [5 £ V £ ]e~ i£ ~ 2ak ~ 1{ - s - z ^y v ^ [E z [6 e Vf\J^ x 0]} (x,x,y)dsdt 
[E z [6 £ V t £ ]e^ 2ak ^ s '~ z ^y^ [E z %V!,\T 2 l e]) {x^y)ds'dA dxdy 
< Ce- 8 J j°° |E {E z [5 £ V t £ }e- i£ ~ 2ak ' lis ~ z)v y^ [E^V^O] } (x,x,y 



dsdt 



E 



[e z [5 £ V £ , 



£] i£- 2a fc- 1 (s'-z)VyV i 



*[E 2 [^^ 2 -^]}(x,x,y) 



ds'dt'dxdy 



+Ce~ 8 J |E[E z [S £ lf]E,[^?]]| |e {e"* ""-^(-^Vy-V* [e^V/]^] ( x>y) 



, e fe-»«H(^)V y .v s [E,^]^ 1 *] (x,x,y)} 



dsdtds dt' dxdy . 



The last inequality follows from the quasi-Gaussian property. Note that in the x integrals above 
the fast variable x is integrated and not treated as independent of x. 
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Let 
and 



g(t) = W 



h{ s ) = e — ~ 2a i°-Hs-* )V y -v a [s £V e T -ig] . 
The same argument as that for Lemma Q yields 

\E[E z [g(tWzMs)}}\ < V 1/2 mg(t)] 2 }^ 1/2 mh( S )} 2 } 

< p( e - 2 (t-z))p( £ - 2 ( S -z))E 1 /2[ 5 2 (t)]]E i/2^2 (s)]) tj8 > z . 

\E[E z [ g (tWz[g(t')]]\ < E^iEMVf^mgit')] 2 } 

< p{E-\t-z))p(e-\t>-z))E 1 %\t)]E l l 2 [ g \t% t,t'>z; 
\E[E z [h(s)]E z [h(s')}}\ < E 1 /2[ Ez [/ i ( s )] 2 ]E 1 /2[E 2 [/ l ( s ')] 2 ] 

< p(e- 2 (s-z))p(e~ 2 ( S / -^))E 1 / 2 [/i 2 (s)]E 1 / 2 [/i 2 ( S / )], s,s'>z. 
Combining the above estimates we get 

/ f-oo \ 4 

?e ^e/i||2 



E||^/:^|^ < 2C / p( S )ds / E[ 5 (z)] 2 E[^)] 2 dxdy 



o 

< 2C ( / p(s)da ) / E[<S e ^/] 2 E 
4 







e -fc- te *- l (-*)VyV a [ 5e y/^-i0] 



dxdy 



o 



< 8C / E[F z e ] 2 / E e~ i£ " ( s ~ z ) v y v * [^F/^-^j (ixdy 



< 8CN p( 8 )<fo) E[V 2 f jE{e^[tf(x,p + q/2)-»(x,p-q/2)] 

x e^" 1 (»-^)P-q/e 2Q t> £ (dq) } 2 dxdp 

< 8C p( a )<k\ E[V Z £ } 2 J [0(x, p + q/2) - 0(x, p - q/2)] 2 $(£, q)d£dxdqdp 

□ 

Eq. I|70|) is convenient for estimating the second moment of p • V x £|jC?.# and £. Z C Z C Z 9 which by 
(|7U|) and (|21() have the following expressions 

^ 2 - 1 {p-V x £^}(x,x,y) 

- 2 V y • v x ^° e - J£ ~ 2Qfe " 1 ( i ~^ v ^ v - {E z [<y 6 lf]e- fc " 2afc " 1 (-*) v ^- Vfi [E^V^^e] } (x,y)dsdt 



IE 



IE 



+IE 



-i£- 2a fc- 1 (t-2)V; 



1 {^^}(x,x } y) 



The same calculation as in Lemma |2] yields the following estimates: 
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Corollary 2. 

E||p.V x £^||2 

< 32C (J p(s)ds^ |E[V y l//] 2 J [V x 0(x, p + q/2) - V x fl(x, p - q/2)] 2 $(£, q^dxdqdp 
+ E[V Z £ } 2 J [V x 0(x, p + q/2) - V x #(x, p - q/2)] 2 |p| 2 <D(£, q)d£dxdqa!p} ; 

E\\£ £ z £ £ z £ £ z e\\ 2 2 < 32C (J^ p(s)ds\ E[F/] 4 ^ [0(x, p + q/2) - 0(x, p - q/2)] 2 $(£, q)d£dxdqdp 

/or some constant C independent of e. 

Let 
(71) 

be the 1-st perturbation of f(z). 
Proposition 1. 



ft(z) = kef'{z)(wiCl9 



lim sup E = 0, lim sup |/f(^)| =0 in probability 



Z<Z 



z<z 



Proof. We have 

(72) 

and 

(73) 



EO/f^O^II/'llcoll^olhEll^lb 



sup |/f(z)| < ell/'HoollWolla sup \\£ e z e\\ 2 . 

Z<Zo Z<ZQ 



The right side of (|72|) is 0(e) while the right side of (|73|) is o(l) in probability by Chebyshev's 
inequality and Assumption 3. 

Proposition ^ now follows from (|72|) and (|73|) . □ 

Set / e (z) = /(z) + ff(z). A straightforward calculation yields 

= e/ , (z)(^,p-V x ^)+6/ ,/ (z)(^,p-V x 0)(iy z e ,^ 
+fc 2 /'(z) <V|,£^\ + fc 2 /"(z) (W* x ,C%0) (W £ z ,£ e z 



-~f{z) (W £ ,£ £ 



l -f'{z) (W £ z ,p- V x 0) + k 2 f'(z) (WI,£ZC%9) + k 2 f"{z) {wic\ 

'wicie 



and, hence 
A £ f(z) 

K, 

(74) +e [/'(*) ^p • V x £^ + f(z) {W £ z ,p ■ V, 

= ^(s) + Ai(z)+Al(z) + R\{z) 
where A\(z) and ^IC- 2 ) are the ^(1) statistical coupling terms. 
Proposition 2. 

lim supE|i?f(z)| 2 = 
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Proof. 



Clearly 



I -fill < e 



>liwbiiiiip- v x fl|| 2 ||^ii2 + ii/'i 



lim supE|i?f(z)| 2 = 0. 



|p-v x (,c* 



£^0 



z<z 



by Lemma [I] and Corollary^ □ 

For the tightness criterion stated in the beginnings of the section, it remains to show 

Proposition 3. {A £ f £ } are uniformly integrable. 

Proof. We show that {^4?}, z = 0,1,2,3 are uniformly integrable. 
For this we have the following estimates: 

\A%{z)\ < i||/lUI|Wo|| 2 ||p-V x 0|| 2 
\A\{z)\ < ^H/'lloollWolbll^^lb 

\A%{z)\ < ^uriuiiwoiiiii^ibii^ib. 

The second moments of the right hand side of the above expressions are uniformly bounded as e — > 
by Lemmas n and 12 and hence Aq(z), A\ (z), A^(z) are uniformly integrable. By Proposition R\ 
is uniformly integrable. □ 

6.2. Identification of the limit. Our strategy is to show directly that in passing to the weak limit 
the limiting process solves the martingale problem with zero quadratic variation. The uniqueness 
of the limiting deterministic problem then identifies the limit. 

For this purpose, we introduce the next perturbations /f , /f . Let 

(75) 

(76) 
where 

(77) Q A o fl)(x.p.y.q) - E 

and 



4 WO = J V(x,p)Qi(^®e)(x,p,y,q)^(y,q)dxdpdydq 
= J Qie(x,p)^(x,p) dxdp, V^GL 2 (M M ) 



£*0(x,p)£*0(y,q) 



Clearly, 
(78) 
Let 

and 

Let 
(79) 

(80) 



q;0(x, p ) = e c e K £ e t e(x,p) 



4%)=e ^,c%e)U,c 



Q 2 (6 0)(x, p, y, q) = E [C £ Z 6( X , p)££0(y, q) 
Q' 2 8(x, P )=e\c £ z C £ z 9(x,p: 



4 2) W = / ^(x,p)Q 2 (0®0) (x,p,y,q)^(y,q) dxdpdydq 
4 2) (V0 = y^(x,p)V(x,p) cfacdp 
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Define 
(81) 

(82) 

Proposition 4. 



fS(z) 



2 1,2 



2 

e 2 k 2 



XVI, C%0) -A%\w e z ) 



f{z) (wicicie)-Af\w, 



lim sup E|/f (*)| = 0, lim sup E|/f (z)| = 0. 



£^0 



z<z 



z<z 



Proof. We have the bounds 



supE|/f(*)| < sup^HnU llWoll^EII^I^+E^f^)] 

Z<Z() Z<ZQ 



sup E | fl (z) | < supe 2 k 



2 1.2 II ,/ 



||^o||2E||^|| 2 +E[4 2) (^)] 



Z<Z0 Z<«0 

The right sides of both tend to zero as e — > by Lemma 1 and 2. 
We have 

^/|(z) = ^/''(^[-(W-^^/^^N+^CWI) + i? 2 (:) 



A £ fi(z) = k 2 f(z) -(w!,cucte))+A^(w : 



with 



2 

-2 ,2 -■// 



(83) 



ww 



^ (w^,p ■ v x (£») + * ( w £ z ,£ z £ z < 
k 



(2) 

where G\ denotes the operator 



Q 2 (6 <g> 9){x, p, y, q)^(y, q) dydq. 



Similarly 



(84) Rl(z) = e 2 k 2 f(z) 



1(W!, P - VACIC16)) + k - ( 117. £:£:£:0 



WW 



wi,c%c%e)-A¥\w*) 



Proposition 5. 



lim sup E|/J|(z)| = 0, lim sup E|U|(js)| = 0. 



e^0 



z<z 



e^0 



z<z 



□ 



Proof. Part of the argument is analogous to that given for Proposition The additional estimates 
that we need to consider are the following. 
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In R £ 2 : First 



sup e 2 E 

Z<ZQ 



e 2 / E [iy 2 £ (x, p)W«(y, q)] E p V x £^(x, p)£|0(y, q) 



dxdydpdq 



< e 2 E[W 2 £ (x,p)W 2 £ (y,q)]E 1 / 2 [p- V x £^] 2 (x,p)E 1 / 2 [£^] 2 (y,q)dxciy(ip(iq 



which is 0(e ) by using Lemma [TJ Corollary [T] and the fact E [Wjr(x, p)Wf(y, q)] G L 2 (M 4d ) in 
conjunction with the same argument as in proof of Lemma 1; Secondly, 



sup eE 

z<z 



sup £ ||Wo||2E||^ 2 V 2 £ || 2 

z<z 



supe||Wo|| 2 E||£ 2 E 

Z<ZQ 



w £ \ 



supellWolbEHJ^- 1 /:^ 

Z<ZQ 



T 2 x l £ z Q®T 2 x l £ z t 



T 2 X W £ Z \\ 2 . 



Let 



We then have 



h 8 = e~ ik l£ 2a ^ s - z ^y^[s £ v z £ ^ 2 h 



= E 
< E 1/2 



t 2 1 c £ z e ® f^ 1 cit 



F 2 l w e zh 



f' f'OO 

J J <5 e V 2 £ (x,y)E [E z [h s {^y)\& z [h t {d^ ,dy')}] F 2 l W £ (x', y')dx! dy' dsdt 
\5 £ V £ ( X ,y)\p(e- 2 (s - z))p(e~ 2 {t - z))¥}' 2 [h s {^y)] 2 



dxdy 



E 1 / 2 [h t (dx',dy')] 2 \J r 2 1 W £ {x',y')\dx'dy'dsdt 



dxdy 



< E 1/2 



.-4 



\5 e V £ (x,y)\p(e- 2 (s - z))p(e- 2 (t - z))E^[h a (x,y)] 2 



E[h t (dx',dy')] 2 dx.'dy' S \ (J |WJ(x', p')|Wdp'J dsdt 



d~xdy 



Recall that \\W £ \\ 2 < \\W \\ 2 and 



J E[ht(dx.', dy')fdx!dy' = f [0(x, p + q/2) - 0(x, p - q/2)] 2 $(£, q)d£dqdxdp 



< oo 



19 



so that 



CP ® £ e z 



< ||Wo|| 2 lE 1/2 ||^iE 



1/2 



2 Wzl|2 

oo 



|5 £ y/(x,y)| /3 (e- 2 ( S - z))p(e- 2 (t - z))E l / 2 [h s (x,y)] 2 dsdt 



ebcdy 



< HWolbis 1 / 2 !!^!!! r^xpiE[5 e "^] 2 ^ ^- 8 ^°° j0 (^- 2 ( s — ^)) j0 ( £ - 2 (t — ^)) 



xp(e~V - z))p(£-V - ^))IE 1 / 2 ||^ s |||M 1 / 2 ||^|||rfs^rfs / (it / 



< ||Ty ||2lE 3/2 ||^||2 fsupE[<5 e V; 
Recall from (|63|) that 



= 12 



OO 



p(s)ds 



< oo. 



E\\h s \\ 2 



Hence 



J [0(x, p + q/2) - 0(x, p - q/2)] 2 $(£, q)d^dqdxdp < oo. 

-0(e). 



sup eE 

Z<Z 



In 



sup eE 

z<z 



which is 0(e) by Corollary [21 
The other two terms in i?| 



^Z ) '^Z^-'Z^z' 



< e||Wo|| 2 supE||^£^|| 2 

z<z 



£ 2 E|<^,p-V x (Q^)>| < e 2 ||^o||2E 1 / 2 ||p-V x E[£^]||J 

< e 2 ||^ ||2||E[p-V x ^£^]|| 2 

which is 0(e 2 ) by Corollary [2 

sM\(W!,LlQ! 2 e)\ < e ||Wo|| 2 E||^E[Z^]|| 2 

< ell^olbEll^-^^E^- 1 ^^]^ 



< e ||Wo|| 2 ( supE 1 / 2 \5 E V*\ 2 ] E 1 / 2 ^^ 2 



which is 0(e) by Lemma |2J 



□ 



Consider the test function f £ (z) = f(z) + f{(z) + /f ( z ) + /f(-z)- We have 
(85) „4 £ / £ (*) 



1 



p • V x 0) + k 2 f"{z)A ( i\wi) + k 2 f'A^(W £ ) + R\(z) + RUz) + R £ 3 (z). 



.2 f/jWi 



Set 

(86) i? £ (z) =^(z)+ J R|(z)+ J R|(z) 

It follows from Propositions 3 and 5 that 

lim supE|i? e (z)| = 0. 
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Proposition 6. 



lim sup sup A^\ip) = 0. 



Proof. We have 



4%) = y^(x,p)Qi(0 0^)(x,p,y,q)V(y,q)dxdpdydq 

= ^ y V(x,p) [Qi(0®0)(y,q,x,p) + Qi(0®0)(x,p,y,q)] V(y, q) dxdp dydq 

where the kernel Q\ (9 <8> 0) can be written as 
Si ® 0) (y, q, x, p) + Qi (9 9) (x, p, y , q) 

/oo /* 
ds / dp / $(s,p / )e ip '-( x - y )/ £2a e-' fe " lsp ' p ' £2 " 2Q [0(x,p + p72) -0(x,p-p72)] 

x [0(y,q + p'/2)-0(y,q-p72)] 
= 2vr y e V-(x-y)/s 2 " ^ (X) p + p / /2) _ e(Xj p _ p / /2) ] 

x [0(y, q + p'/2) - 0(y, q - p'/2)] ^fc^p • p'e 2 - 20 , p')c*p' 



which is uniformly compactly supported on R . For smooth and compactly supported <3>, Qi(9 <g> 
0) tends to zero fast than any power of e uniformly outside any neighborhood of x = y while 
stays uniformly bounded everywhere. Therefore the L 2 -norm of Q\{9 <8> 9) tends to zero and the 
proposition follows. □ 



Similar calculation leads to the following expression: For any real-valued, L 2 -weakly convergent 
sequence if) £ — > ip, we have 

lim A?\i/>) 

POO P 

= lim / ds / dwdqdxdp ip e (jc,p)^(w,q)e isw e~ ik ' lsp - Ci£2 ' 2a 
£ -°Jo J 

e -ifc-l S |q| V-»-/2 p( X) p + q ) _ 0( x> p )] _ e ik-^\'e'-^/2 ^ p) _ ^ p _ q) 

= lim [ ds [ (iqdxdp V e (x,p)6(s,q)e^ ife_lsp ' qe2_2a 
Jo 7 

-ifc-^lql V" 2 «/2 [0 (Xj p + q ) _ 0( Xj p )] - gifc-^NI V- a -/2 [e(X) p) _ e(x> p _ q) 



Note that the integrand is invariant under the change of variables: 



s^-s, q^-q. 
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Thus we can write 
lim^WO 

£—♦0 



lim — 

e^o 2 



<h I c/qr/xf/p r ; (x.p )<]>(. s-.qje - ** ls P^ 2 " 2a 

ifc- 1 s|q| 2 e 2 - 2o! /2 



CO 

e -ifc-i 8 |q|V->-/2 ^ (X) p + q ) _ 0( X; p) ] _ e * 



[0(x,p)-0(x,p-q)] 



limvr y dqcixdp V e (x,p){$(e 2 ~ 2Q fc- 1 (p + q/2) -q,q) [0(x, p + q) - 0(x, p) 

-^(^-^fc-Hp - q/2) • q,q) [0(x,p) - 0(x,p - q)]} 
liniTr / dqdxdp^(x,p)(cD( £ 2 - 2 " |q| ~ |p| ,q-p) [0(x, q) - 0(x, p)] 



£^0 



2A* 



-^(£ 2 ' 2 jp|2 2 J q|2 ,P " q) [<?(x, p) - 0(x, q)] | 



Clearly, we have that 



I |2 I 12 

lim 2,7 / dqdxdp ^ £ (x,p)$(£ 2 - 2a lPI 2 ~' ql ,p- q) [0(x,q) - 0(x,p)] 



lim / dq[0(x,q)-0(x,p)]$( £ 2 - 2 «i q ^— ^,q-p) 
£— »o J zk 

$((2fc)- 1 (|q| 2 -|p| 2 ),q-p), a=l 
dq[0(x,q) -0(x,p)] X { $(0,q-p), 0<a<l 

3, a > 1 



in L 2 (E 2d ). Therefore 



lim A\ '(tp 6 



$((2A;)- 1 (|q| 2 -|p| 2 ),q-p), a=l 



dqdxdp^(x,p)[6(x,q) -0(x,p)] X { $(0,q-p), 

0, 



< a < 1 
a > 1 



Am 



Recall that 

(87) M|(0) = f(z)- I .4 : /--(.s)r/.s 



f(z) + ff(z) + f!(z) + f!(z) - J* ~f'(z) <WJ,p ■ V x 0) ds 



f"(s)4\w!) + f'(s)A?(W! 



ds 



R £ (s)ds 



is a martingale. The martingale property implies that for any finite sequence < Z\ < Z2 < z^ < 
■ ■■ < z n < z, C 2 -function / and bounded continuous function h with compact support, we have 

(88) E{h «w£ , 0> , (w* , 0) , (w! n ,e)) [m s z+ m - Mm] } 

Vs > 0, Z\ < Z2 < • • • < z n < z. 

Let 



0, 



Af(z) = f'(s) 



k 



(W z ,p- V x 0) + A: 2 ii(W 2 ) 



In view of the results of Propositions 11121 151 IH 151 we see that f 6 (z) and A 6 f £ (z) in (fHTj) can be 

replaced by f(z) and Af(z), respectively, modulo an error that vanishes as e — > 0. With this and 
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the tightness of {WJ} we can pass to the limit e — > in (|88j). We see that the limiting process 
satisfies the martingale property that 



E{h((W Zl ,9) , (W Z2 ,e),...,(W Zn ,0)) [M z+S (6) - M z {0))} = 0, Vs > 0. 



where 



M z (e) = f(z)- / Af{8) da. 
Jo 

Then it follows that 

E[M Z+S (9) - M z (6)\W u ,u < z] =0, Vz,s>0 
which proves that M z {6) is a martingale given by 



M z {6) = f(z) - / f'(s) 



(90) 



Choose f(r) = r and r in (J9(J|) we see that 



^ <W s ,p- V x 0) + fc 2 Ai(W s 



(is. 



AfW(e) = <Wi,6>)- 

is a martingale with the null quadratic variation 

M«(0),M«(0) 

Thus 



/(*) " / /'(«) 



<W s ,p- V x 0) + A; 2 Ai(W fl 



d s = /(0), Vz > 0. 



Since (WJ,#) is uniformly bounded 

\(w!,e)\<\\Woh\\o\\ 2 

we have the convergence of the second moment 

YimE{(W z s ,e) 2 } = (W z ,e) 2 
and hence the convergence in probability. 

7. Proof of Theorem 2 
7.1. Tightness. Instead of Q6UJI we use the following the corrector 
(9l4^(x,x,p) 

| w c-w.'v-Xx, p + ^ _ „ (x , p _ ."wpSJw 

which satisfies the corrector equation (|65l) . Its inverse Fourier transform is given by 

e _, £ - 2Qfc -i (s _, )Vy .V x [ Ez (X)y)ds . 

Instead of Lemma ^ Corollary ^ Lemma |2] and Corollary |2] we have 
Lemma 3. Assumption 1 implies that 

roo 1 2 p 

[q-V p 0(x,p)] 2 <D(C,q)d^q 

which has a bounded support and is bounded uniformly in x, p. 
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(93) 


lim sup E 

e^0 


LP 


2 (x,p) < 


/ p(s)ds 
Jo 


7 



Corollary 3. 



£^0 













2 (x,p) < 


/ p(s)ds 
Jo 


J [P'V x [q 



[p-V x [q- Vp0(x,p)]]^(£,qRdq 



which has a bounded support and is bounded uniformly in x, p. 
Lemma 4. 



limsupE||£^0||! < 8C 

limsu P E||£ £ £ £ 0||| < 8C 

/or some constant C. 
Corollary 4. 

limsupE||£*£*£*0||| < 32C 

£->0 

limsupE||p • V^C £ z C £ z 9\\l < 32C 

e^O 



p(s)ds nVzY / [q • V p 0(x ) p)]^(^,q)dedxdqdp 



p{s)ds E[V Z ] 2 / [q • V p 0(x,p)] 2 $(£,q)d£rfxdqdp 



p(s)ds E[V z f / [q • V p 0(x, p)] 2 $(C, q)dedxdqdp; 



|E[V y K] 2 y [q • V p V x fl(x, p)] 2 $(£, q)d£dxdqdp 
+ E[y 2 ] 2 y [V x q • V p 0(x, p)] 2 |p| 2 <l>(£, q)d£rfxdqdp 



for some constant C . 

The rest of the argument for tightness proceeds without changes. 

7.2. Identification of the limit. With straightforward modification on the estimates of the re- 
mainder terms -Rf, an d R% the same argument for passing to the limit e — ► as before carries over 
here. 

In particular, Proposition |H1 can be proved as follows. 



Proposition 7. 



lim sup A^(ip) = 0. 

e^O iu/,11 i 



Proof. The kernel Q\{9 <8> 9) has the following expressions 
Gi(0®0)(x,p J y J q) 

roc p 

= Hs, p')e tp '-( x -y)/ e2Q e- ifc " ls P-P' e2 " 2Q e 2a - 2 [0(x, p + e 2 - 2 V/2) - 0(x, p - £ 2 - 2 V/2)] 

xe 2a ' 2 [9(y, q + e 2 ~ 2 V/2) - 9(y, q - e 2 - 2 V/2)] dp' ds 

= lim IX f $(A; _1 p • p'^ 2 ", p ') e iP'-(x-y)/^ £ 2a-2 p + £ 2-2a p , /2) _ (X) p _ £ 2-2a p / /2) l 

xe 2 °- 2 [0(y, q + e 2 - 2a p'/2) - 9(y, q - e 2 - 2 V/2)] dp'. 
Note that 

e 2a " 2 [0(x, p + e 2 - 2a P 72) - 0(x, p - e 2 - 2 V/2)] — p' • V p 0(x, p) in C™(R 2d ). 



Thus the L 2 -norm of Q\( 
tionEl 



tends to zero for the same reason as given in the proof of Proposi- 

□ 
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To identify the limit we have the following straightforward calculation: For any real-valued, 
L 2 -weakly convergent sequence ip £ — > ifi 

lim A { ^{^ £ ) 



= lim/ ds dwdqdxdp 4; £ (x,p)<5>(w,q)e isw e- ik ~ lspci£2 ~ 2a e 4a -* 
J 

e -ifc-l a |q| V-*-/2 [0 (X) p + £ 2~ 2aq) _ p) ] _ e <fc-l S |q| V-*«/2 p) _ 0(x> p _ £ 2-2a q) ] ' 

roc /* 

= lim / (is / dqdxdp ^ e (x,p)6(s,q)e- i ^ lsp q£2_2Q e 4Q - 4 
Jo J 

e- 4fc -Mq|V--/2 ^( X) p + e 2-2a q) _^ (X)p) ] _ e ^ S |q|V- 4 "/2 [0( X) p) - 0(x, p - e^q)]" 
= ^ y dqdxdp i/>(x, p)*(0, q)(q ■ V p ) 2 0(x, p) 

= ^iW 

for a £ (0, 1). For a = 1 we have the same result as in Theorem 1 Case (ii); for a > 1, the limit is 
identically zero. 

8. Proof of Theorem 3 

The proof of the result for Case (i) and (iii) is identical to that for Theorem 1, Case (i) and (hi), 
respectively. So in the sequel we focus on the second case a > (3. 
Introducing a new parameter 

6 = 6? 

we can rewrite the equation as 

f)W £ r> k 

(94) ^ + £ • VWJ + -C £ Z W £ = 
with 

(95) £ e 2 T^(x,x,p) 

= ii 1 ~ a ^ J e iq '* [WJ(x, p + q/2) - W 2 £ (x, p - q/2)] V £ (dq), x = X6~ 2a ^ , /3 < 1 

with 

(96) V z £ (dq)=V(^,dq), V £ {^) = V{^). 
Note again that 

(97) T^C%e = -ifW^^y^fl 
with 

(98) 5 e t?(x, y) = Kf(x + y/2) - F/(x - y/2). 

We will work with (|95|) and (|96|) and construct the perturbed test function in the power of e. 
First we note that 

(99) E [C £ z 0] 2 (x, p) = 6 2 ~ 2a ^ j [0(x, p + q/2) - 0(x, p - q/2)] 2 <D(£, q)d£dq 

which has an e-uniformly bounded support. 
Instead of (|6(Jj) we define 

£^(x,x,p) = ¥ ^ y^e^*"^-')^^ 
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with 



x = xe 



-2a//3 



After the partial inverse Fourier transform 



n+a/0 



(100) J7 1 £^(x,x,y) 

£ ' ■ j z 
The corrector equation holds again: 



(101) 



[0(x, p + q/2) - 0(x, p - q/2)]' q)d£dq 



Following the same argument as in the proof of Theorem 1 we have the following estimates: 

Lemma 5. 

r 1 2 r f°° 

(10*1 £*0 (x,p) < e 2 ~ 2a/f3 / p(s)ds 

which has an e-uniformly bounded support. 
Corollary 5. 

(103) E p • V x £^6> 2 (x,p) 

r r°° "i 2 /• 

< J Q P( s ) ds J [p-V x 0(x,p + q/2)-p-V x 0(x,p-q/2)] 2 $(e,q)^q 

which has an e-uniformly bounded support. 
Lemma 6. 

E||£|£>||| < &-*»lf*%C ( [°° p(s)ds) E[V Z } 2 [ [0(x, p + q/2) - 0(x, p - q/2)] 2 <l>(£, q)d£dxdqdp 



E\\£&\\ 2 2 < e 4 - 4Q//3 8C p( s )tfe^ E[V Z } 2 J [9(x, p + q/2) - 0(x, p - q/2)] 2 <£(£, q)d£dxdqdp 



/or some constant C independent of e. 
Corollary 6. 



E||p-V x £^|| 2 

< e A ~ Aa ^32C (J °° p(a)d^ |E[V y y 2 ] 2 y [V x 0(x, p + q/2) - V x 0(x, p - q/2)] 2 $(£, qXdxdqdp 
+ E[V Z } 2 J [V x fl(x,p + q/2) - V x 0(x,p-q/2)] 2 |p| 2 <lK£,qRdxdqdp} ; 



< e 6 - 6a/l3 32C p(s)ds^J E[V Z ] A J [0(x, p + q/2) - 0(x, p - q/2)] 2 $(£, q)d£dxdqdp 

for some constant C independent of e. 

As a consequence of the divergent factor e l ~ 1 /^ in the above estimates the previous proof of 
uniform integrability of A e [f(z) + iff] (e.g. Proposition |3J) breaks down. For both the tightness 
and the identification we shall use the test function 

f%z) = f(z) + ft(z) + fl(z) + fl(z) 
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where 
(104) 

(105) 
(106) 



/!(*) 

/!(*) 



kef f (z)(W!,C £ z 



;2 



A: 2 



/'(*) (wicicie)-Af\w, 



(2) (2) 

with \A 2 > given as before. 

Following the same procedure as in the proof of Theorem 1 we obtain 

(107)A £ / £ (^) 



\f\z) {Wl,p- V x 9) + k 2 f"{z)A ( i\wl) + k 2 f'A?\wl) + Rl(z) + B&z) + A%{z). 



where 

(108) R £ 2 {z) 



2 

^2il2 /■/// 



+e^f"{z){w* z ,c%e 

-e 2 k 2 f"(z) 



(109) R%{z) = e 2 k 2 f\z) 



|(n':. 



.fc 2 
2 



A {2) {Wt) 



e 2 k 2 f'(z) 



and A%, A [ 2 l) , A*p , G {2) , Q' 2 all have the same expressions as in the proof of Theorem 1. 

With the assumption (3 > 3/4, Propositions and |S] hold true. Let us remark that the most 

severe terms due to the divergent factor are 



(110) 
(111) 



sup eE 

z<zo 

sup £E 

z<z 



Wf,£lG {2) Wf 



Wl£l£%£%t 



(cf. Corollary EJ) . 

To satisfy we need to show 



Proposition 8. 



lim sup \ ff x | =0, j = 2,3, in probability. 



z<z 



Proof. We have the estimates 

sup|/f )Z | < supe 2 ^ 2 ||r||oo[||Wo|||||^||| + 4 2) (^)" 

z<zo z<zo L 

sup \flj < supe 2 A ; 2 ||/'|| 00 [||Wo||2||^^||2 + 4 2) (^ 

z<zo z<zo L 

which vanish in probability by using Assumption 3, Lemma El El and Chebyshev's inequality. 
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□ 



Proposition 9. 



lim sup sup A^\ip) = 0. 



£^0 



Z<ZQ \\lh\\ 2 = l 



Proof. The kernel Q\{9 <g> 9) has the following expressions in terms of e (not e): 



Qi(0®0)(x,p,y,q) 

= e 2 ^ 2a J°° J $(s, p >) e ip'i*-y)e- ik - ls ve'c 2f3 - 2a [0(x, p + p'/2) - 0(x, p - p'/2)] 

x[e(y,q + p72)-(9(y,q-p72)] dp' ds 
= s 2f3 - 2a 7r J ^(k _1 p • p'e 2/3 - 2Q , p')e ip, -( s -^ [0(x, p + p'/2) - 0(x, p - p72)] 
x [<9(y,q + p72)-e(y,q-p72)] dp'. 



Taking the L 2 -norm and passing to the limit we have 



lim / dxdpdydci\Qi(9 tg> 6)(x, p,y,q) 

e->0j 

lim 7r 2 J dxdpdydq J 6( — ^ 



$(w,p')dw 



ip'-(x-y)/ £ 2 



[0(x, p + p72) - 0(x, p - p72)] [0(y, q + P 72) - 0(y, q - P 72)] dp 



/|2 



lim7r / dxdptiydq 



y y~j J &(w,p±)dw 



3 ipx-(x-y)/e 2 



X [0(x, p + p ± /2) - 0(x, p - p±/2)} [0(y, q + p ± /2) - 0(y, q - p ± /2)] dpj 



where pj_ • p = 0, pj_ G M rf 1 . In passing to the limit the only problem is at the point p = 0. But 
the integrand in the above integral is bounded by 



(112) 



C|p| 2 , some constant C < oo 



which is integrable in a neighborhood of zero if d > 3. Hence the L 2 -norm of Q\ (9 ® 9) tends to 
zero by the dominated convergence theorem. □ 
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We have the following straightforward calculation: For any real-valued, L 2 -weakly convergent 
sequence ip £ — > ip, 

lim a[ 1] W) 



f'OO r 

lim / ds I dqdxdp ip e (x, p)l>(s, q)e 2/3_2a 
Jo J 

e -ifc-l s (p +q /2). qe 2/3-2a ^ X)p + q ) _ 0( Xj p)] _ e -ife- 1 Kp-q/2)-q£ 2 ' 3 - 2a [0( X) p) - 0(x, P - q)] 

lim vr / dqdxdp V £ (x, p) [*(Jfc -1 (p + q/2) • qe 2/3 ~ 2a , q)e 2/3 " 2a [0(x, p + q) - 0(x, p)] 
e^o y L 

-*(Af X (p - q/2) • q^-^.q)^- 20 [0(x,p) - 0(x,p - q)f 



lim 7T / <iq<ix<ip -^ £ (x, p) 



I |2 I 12 

IgJ -IPI £ 2/?-2a >q _ p ) e 2 / 3-2a p + q ) _ g( X) p ) 



2A' 



|q| 2 ~ IpP ^20-2a 



2k 



£ 2 ?- 2a , q - p)e 2 ^- 2Q [0(x, p) - 0(x, p - q)] 



2tt / dqdxdp ^(x, p)<5( 



|q| 2 -|p| 2 , 
2k ' 

= 

following from the strong convergence of 



J q — p)dw 



[0(x,p + q)-0(x,p) 



/ ^ ^ ~ P) £2/3 ~ 2 ° ^ P + ~ °^ P)l 

-^( |q|2 2 ~ |p|2 g 2/3 ~ 2a , q - p)^- 2Q [0(x, p) - 0(x, p - q)] 



dq, V$ E C, 



o°(R 2d ) 



in L' 



Proposition 10. A £ f e (z) is uniformly integrable. 

This, of course, follows from the fact that each term in Q1U7JI has a uniformly bounded second 
moment. Therefore we have completed the tightness argument. Moreover, we have also identified 
the limiting equation. 



Introducing a new parameter 
the fast variable 



9. Proof of Theorem 4 



e = e h 



x = xe 



-2a 



and the rescaled process 
(113) 

we rewrite the equation as 
(114) 



y/(dq) = V(^,dq), y/(x) = V(±x) 



m * + ? • vw! + W = o 



dz k 
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with, in the case (i), 

(115) £^ 2 £ (x,x,p) 

= i J e^V a - 2 [^(x,p + e 2 - 2Q q/2)-^(x,p- £ 2 - 2a q/2)] V*(dq), 

and in the case (ii) 

(116) £ e 2 ^(x,x,p) 

= ie l ~ a l p J e iq 'V Q - 2 [W 2 e (x, p + e 2 - 2Q q/2) - W*(x, p - e 2 - 2o q/2)] V; £ (dq). 

Taking the partial inverse Fourier transform we get in the case (i) 

(117) ^-i£|0( X) x )y ) = -i ( 5 £ y/(x,y)^ 2 - 1 0(x,y) 
and in the case (ii) 

(118) jr-l£^( x , X ,y) = -ie 1 -^,5 £ V; £ (x,y)^ 2 - 1 0(x,y) 
with 

(119) W(x,y) =e 2a ~ 2 [F/(x + e 2 - 2a y/2) - V; £ (x-e 2 - 2Q y/2] . 

The proof for the case (i) is entirely analogous to that for Theorem 2 and we will focus on the 
case (ii) in the sequel. And we will work with (|116|) and (|118|) and construct the perturbed test 
function in the power of e. 

First we note that 

(120) hmsupe- 2+2a ^E[/:^] 2 (x,p) = /[q-V p 0(x,p)] 2 <i>(£,qKdq 

which has an e-uniformly bounded support. 

As in the proof of Theorem 3, we carry the analysis in the power of e = e@ . We consider the 
rescaled process (|9*6^l and its sigma algebras. 

We introduce the corrector 



£p>(x,x,p) 

i 



j°° j e iq.x e ifc- 1 ( S -,)p.q/ £ - 2 -/' 3 e 2a-2^ (X)p + £ 2-2a q/2) _ ^ p _ £ 2-2<* q / 2 )]E|y/(dq) 



which after the partial Fourier inversion becomes 



(121) ^£^(x,x,y 



=l+a//3 



£ 



j ifc- 1 ( S -«)VyV x /£ 2 «// 3 w err T/ ei T -l 



The corrector solves the corrector equation eq. ()1U1|) . 

Following the same argument as in the proof of Theorem 1 we have the following estimates: 



Lemma 7. 

(122) limsup;r 2+2a//3 E \c £ z 

which has a compact support. 
Corollary 7. 



(x,p) < 



p(s)ds 



[q-V p 0(x,p)] 2 <I>(£,qRdq 



limsupe" 2+2a//3 E 
which has a compact support. 







" />oo 




P • V x £ e ^ 


2 (x,p) < 


/ P(s)ds 
Jo 
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[p- V x q- V p 0(x,p)] 2 $(£,q)d£dq 



:-!() 



Lemma 8. 

limsupe~ 4+4a//3 E||^£^||^ < 8C (J p(s)ds\ E[V Z } 2 J [q • V p #(x, p)] 2 $(£, q)dfdxdqdp 
limsupe- 4+4a//3 E||£f£ £ Jll9 < 8C 



LSUpf 

£^0 



oo \ 4 

— |I2 _ - !/ p(s)<fej E[K] 2 /[q.Vp0(x,p)f*(Cq) f 7^/xJcK/p 



/or some constant C independent of e. 
Corollary 8. 

limsupe- 4+4a//3 E||p • V X £*£*0||1 



< 32C (J p{s)ds^j |E[V y F z ] 2 y [V x q • V p #(x,p)] 2 $(£,q)ci£dx(iqc2p 
+ E[K] 2 y [V x q • V p #(x, p)] 2 |p| 2 $(£, q)d£dxdqdp 



lim sup £-6+6a//3 E 1 1 c E z L e z C%e 1 11 



< 32C< (y p(s)ds^ E[y 2 ] 4 y"[q- Vp^(x,p)] 2 $(^q)(iedxdq(ip 

/or some constant C independent of e. 

The rest of the argument follows the general outline of that of Theorem 3 Case (ii). 
Let us know verify that the quadratic variation vanishes in the limit. 



Proposition 11. 



lim sup sup A2 (-0) = 0. 



>z<z \\ip\\ 2 =l 

Proof. The kernel Q\{9 <S> 0) can be calculated as follows. 
gi(0®0)(x,p 5 y,q) 

_ _2/3-2q 



xe 2a ~ 2 [6(y, q + e 2 - 2 V/2) - 9(y, q - £ 2 - 2 V/2)] dp' 

7T y £ 2/3_2a $(A; _1 p • p'e 2 ' 3 - 2 ", p ') e V-(x-y)/^ £ 2a-2 p + £ 2-2a p / /2) _ ^ p _ £ 2-2cy /2) ] 



2a-2 



[9(y, q + e 2 - 2a P 72) - 9(y, q - e 2 - 2 V/2)] dp' 



whose L 2 -norm has the following limit 



lim 7r 



&lpl 



p ± p=0 



$(w, p±)dw 



= ipj_-(x-y)e 



|q- V p fl(x,p)| 2 dp, 



dxdp = 



if d > 3 by the dominated convergence theorem because the integrand is bounded by the integrable 
function H112|) in a neighborhood of zero. □ 
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To identify the limit, we have the following straightforward calculation: For any real-valued, 
L 2 -weakly convergent sequence ip 6 — > ifi, 

lim A^\t/) e ) 

= lim r ds ( dqdxdp V £ (x,p)6(s,q)e- i ^ lsp - qe2 ' 3 ~ 2Q e 2/3 - 2a e 4a - 4 
Jo J 

e -ik-^e^^V2 ^ (Xj p + £ 2-2a q) _ ^ p) ] 



(x,p-e 2 - 2Q q) 



lim / ds I dqdxdp V £ (x,p)6(s,q)e- ife ^P-q^V^V " 4 



e -ik-^e^-^/2 ^ (Xj p + £ 2-2 Qq) _ 0(X) p) ] 
_ eifc - ls | q |2 £2+2 ,-4 Q/2 ^ (Xjp) _ 0(X)P _ £ 2-2 aq) ]- 

*(fc _1 (p + e 2 - 2 °q/2) • qe 2/3 - 2Q , q)e 2 - 2Q [0(x, p + £ 2 ~ 2a q) - 0(x, p)] 
$(fc _1 (P - e 2 - 2Q q/2) • qe 2/3 - 2Q ,q)e 2 - 2Q [0(x,p) - 0(x,p - £ 2 ~ 2a q)]l ^ e (x, p)dqdxdp 



lim tt / e 2a - 2 e 2f3 - 2a 



vr^q- V p (JCfc^p-q) / 



TT / fc|p| 1 



• q) 



q- v p 



0(x, p)V ; (x, p)dq<ixdp 



(P± • V p ) 2 6'(x,p)V'(x,p)£ip_ L dx(ip 



where pj_ G 



pd-1 



p_L • p = 0. 



Case (iii): a = (3. 

lim 

e-+0 

= limvr J dqdxdp V £ (x,p)e 4Q " 4 [$(Jfe -1 (p + e 2 ~ 2o q/2) q,q) [0(x, p + e 2 ~ 2a q) -0(x,p)] 

- <f>(k-\p - E 2 - 2a q/2) ■ q, q) [0(x, p) - 0(x, p - £ 2 - 2Q q)] ] 
= tt J dqdxdp ip(x, p)q • V p [^(k^p ■ q, q)q • V p ] 0(x, p) 
= A^). 
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